A smooth cascade of wrinkles at the edge of a floating elastic film 
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The mechanism by which a patterned state ac- 
commodates the breaking of translational symme- 
try by a phase boundary or a sample wall has been 
addressed in the context of Landau branching in 
type-I superconductors [T], refinement of mag- 
netic domains [2], and compressed elastic sheets 
We explore this issue by studying an ultra- 
thin polymer sheet floating on the surface of a 
fluid, decorated with a pattern of parallel wrin- 
kles. At the edge of the sheet, this corrugated 
proflle meets the fluid meniscus. Rather than 
branching of wrinkles into generations of ever- 
smaller sharp folds [3], we discover a smooth cas- 
cade in which the coarse pattern in the bulk is 
matched to fine structure at the edge by the con- 
tinuous introduction of discrete, higher wavenum- 
ber Fourier modes. The observed multiscale mor- 
phology is controlled by a dimensionless param- 
eter that quantifies the relative strength of the 
edge forces and the rigidity of the bulk pattern. 

When a thin rectangular sheet floating on the surface 
of a pool of hquid is compressed along two opposing 
edges, it forms a pattern of wrinkles parallel to these 
edges. Unlike the Euler buckling of an unsupported piece 
of a paper, where the largest possible wavelength is se- 
lected, the wrinkles form at a wavelength A ^ W, the 
width of the rectangle in the direction of the compression. 
Two principles are essential to understanding the ampli- 
tude and wavelength of this pattern: first, a thin sheet 
can, to a first approximation, be treated as inextensible, 
so that the length of a line in the compression direction is 
preserved. Consequently, the wavelength and amplitude 
of the wrinkles are related. Second, the bending energy 
of the sheet favours long wavelengths (and therefore large 
amplitudes) whereas the gravitational energy of the liq- 
uid subphase favours small amplitudes (and hence small 
wavelengths). Thus, the wavelength is selected by a com- 
promise |5| between the bending energy of the sheet and 
gravitational energy (or more generally, any mechanism 
by which the subphase resists distortion). 

We show in Fig[T] a controlled experimental realisa- 
tion of this situation with a polystyrene (PS) sheet of 
dimensions W x L = 3cm x 2cm, and thickness t float- 
ing on the surface of water, and initially experiencing an 
isotropic tension « 7, the liquid-vapour surface tension. 
As the sheet is compressed by a distance A (such that 
A = A/W ^ 1), parallel wrinkles develop in the bulk. 
The one-dimcnsional pattern of wrinkles in the bulk of 
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FIG. 1: (A) Image of a wrinkled PS sheet floating on the sur- 
face of water, compressed along its length between two razor 
blades. (B) Sketch of geometry. (C) Bulk wavelength of wrin- 
kles, A = 27r/go as a function of film thickness, t. The solid 
line is a fit to f"^''*, showing agreement with the prediction of 

the sheet, characterised by a height field C(?/)i has for 
small amplitudes, an energy per unit surface area: 
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The first two terms in u{y) represent bending energy of 
the sheet and the gravitational energy of the fluid, re- 
spectively, while the third term enforces the constraint 
of inextensibility in the limit of small amplitudes, with 
the Lagrange multiplier a being the stress ayy that must 
be applied in the y-direction at the compressed edges, p 
is the density of the fluid, and the bending modulus is 
denoted by S = Et^/{12{1 ~ A^)) where E is the 
Young's modulus and A is the Poisson ratio. The sur- 
face tension does not appear in the energy functional 
shown above due to the fact that the bulk pattern has 
translational symmetry in the ^-direction. However, it is 
important to note that the sheet still experiences a ten- 
sion a-xx ~ 7. Minimising this energy density leads to a 
pattern 



ay) = -VAsmiqy) (sl); ayy = {Bq'+pg/q') (b)(2) 
where the wavenumber q = qo = {pg/ BY^^. With q = qo, 
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FIG. 2: (A) Wavenumber q{x) as a function of distance x 
from the edge of the sheet. The lengthscale over which the 
decay occurs does not change strongly with thickness t of 
the sheet. Exponential fits (solid lines) to q{x) show sys- 
tematic deviations but yield a decay length of 1.8 ± 0.2mm. 
By comparison, the capillary length Ic = 2.7mm (B) Scaled 
wavenumber q{x)/qo versus the scaled distance from the edge 
xqo. Data collapse is good at small and large xqo, but there is 
poor collapse at intermediate distances indicating a multiscale 
evolution. The solid line shows an exponential fit. 



Uyy — —2^Bpg. As shown in FiglH this correctly de- 
scribes the scaling of the wavelength of the wrinkles in 
the bulk. This scaling 5] has been experimentally tested 
0, iS] , and more broadly, applied in situations where the 
wavelength is determined by balancing the bending en- 
ergy against substrate elasticity i^, capillary forces [lo| . 
and stretching under tensile forces 0. 

However, in addition to the pattern in the bulk, in- 
spection of FiglT] shows a striking phenomenon: the par- 
allel wrinkles of the bulk with wavenumber qo, give way 
to a much finer structure of wrinkles near the uncom- 
pressed edge of the sheet. It is this cascade to ever-higher 
wavenumber that we examine in this Letter. 

In Fig. [21^, we show for PS films of thickness t rang- 
ing from 85 to 246nTO, the increase in wavenumber q(x) 
as a function of distance x from the uncompressed edge 
of the sheet. Since the bending modulus B ^ t^, this 
represents a broad range of B. For all thicknesses, as 
the edge is approached, the wavenumber increases to a 
value at the edge qe, that is 2-5 larger than the bulk value 
qo- The evolution to higher wavenumbers occurs over ap- 
proximately the same distance from the edge: though an 
exponential fit shows systematic deviations, such a fit es- 
timates the penetration length of the edge into the bulk 
to be 1.8 ± 0.2mm. 
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FIG. 3: (A) The energy per length of the meniscus, Ucap{qe) = 
2^{K/q^)yJ {pg/^) + qi- (B) Energy per area of a sheet with 
wavenumber q, u{q) = 2{A/q^)B{q'^ + qt). 

Intuitively, a higher wavenumber at the edge is to be 
expected. The fluid meniscus follows the contour of the 
edge of the sheet. To minimise the surface energy of the 
air-water interface it is therefore favourable to reduce the 
amplitude of the wrinkles at the edge. In order to achieve 
this while preserving inextensibility, the wavenumber in- 
creases. This cascade to finer wrinkles must terminate 
at a finite wavenumber, ge, where the gain in surface en- 
ergy is offset by the increased energetic cost of bending. 
Notwithstanding these plausible arguments, previous ex- 
periments in this geometry Q did not show a marked 
effect at the boundary of a film, and find that qe ~ qo- 
We thus need to address some obvious issues: What are 
the relevant parameters that dictate whether a cascade 
should be anticipated? When a cascade is observed, what 
governs the amplification of the bulk wavelength, and the 
length over which the cascade occurs? 

In order to understand why wavenumber amplifica- 
tion occurs in our experiment, we estimate the en- 
ergy cost of a wavenumber q^ at the edge. The cap- 
illary energy associated with the meniscus is Ucap — 
2^{K/ql)yJ{pg/^) + ql, which as sketched in Fig. [4j\, 
is a decreasing function of qe. To compare Ucap, which 
is an energy per unit length, to the energy cost per unit 
area of the affected part of the sheet (Fig. |4j3), we re- 
quire an understanding of the lengthscale over which the 
pattern at the edge penetrates into the sheet and matches 
the bulk pattern. 

This leads us to consider the effect on the energy of 
the sheet of imposing a wavenumber qe at its edge. The 
effect of broken translational symmetry in the ai-direction 
is incorporated by modifying Eq. [l]to: 
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-2A] . (3) 



The new, third, term, ut quantifies the energy cost of 
the deformations in the x-direction under the tension 
7 of the fluid tugging at the uncompressed edge of the 
sheet. When the compressive force in the y-direction is 
much smaller than the tensile force in the x-direction, 
e = 7/cr = 7/V pgB <g; 1, the minimal energy profile is 
determined by balancing bending and tensile forces [see 
Supplementary Information]. We thus develop a length 
scale, Zp, that will govern gradients in the x-direction. 
From Bq^ ^ 7C^j we obtain Ip « pg, which is the 
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capillary length, Ic- This is consistent both with the mag- 
nitude of the typical penetration length found in Fig. [2l 
and with the insensitivity of this length scale to the thick- 
ness t. Furthermore, the stress ratio e ranges between 
6 X 10"'* and 3 x 10"'^ in our experiments, thus validat- 
ing the regime assumed in the foregoing argument. 

A self-consistent argument for the strength of the wave- 
length amplification may be constructed by assuming 
that a single penetration length, Ip, governs the pene- 
tration of the edge wavenumber into the bulk; an esti- 
mate of the bending energy cost near the edge of the 
sheet is then obtained as Uedge ^ lpu{qe) ^ Bq^{A/q'^). 
Setting the edge and meniscus energies to be compa- 
rable, Uedge ~ Ucap We obtain Qe ^ (e)^^^^ {-J / B). 

This argument shows that for e <C 1 (as in our exper- 
iments) the wavenumber amplification is a large, non- 
perturbative effect on the bulk pattern. However, the 



estimate Ue. 



derived above is wrong as it re- 



lies on the oversimplified picture of a single penetration 
length, unlike the cascade observed in Fig. 2 which in- 
dicates an elastic instability that leads to a symmetry- 
breaking sequence of period fissioning into intermediate 
wavenumbers (and correspondingly a sequence of pene- 
tration lengths). A full solution of the nonlinear problem 
llj reveals that this instability results from a logarithmic 
dependence of Uedge log qe , which lowers the energetic 
cost for sufficiently large values of qe- 

Assuming a pattern that is the superposition only of 
two wavenumbers qo and qe , the penetration depth Ip ex- 
hibits a strong dependence on e, ranging from Ic for e <C 1 
to q^^ for e 3> 1 (see Supplementary material). The lat- 
ter limit, together with an energetic estimate similar to 
the one presented above, shows that for e 3> 1 the edge 
effect is only a small perturbation to the bulk pattern 
Q- For e ^ 1 when there are wavenumbers q interme- 
diate between qoand qe, one can show that the length 
lp{q) of the transition zone between wavenumbers q and 
qe, strongly depends on q, ranging from Ic for q ^ qo to 
q^^ for q ^ qe- This is because a periodic wrinkling pat- 
tern with wave number q > qo requires enhancement of 
the compression (t{x) towards the edge (thus increasing 
e close to the edge). That the description in terms of a 
single penetration length Ip is simplified can be seen from 
Fig. 2B, where we show the scaled wavenumber q(x)/qo 
vs. the scaled distance from the edge xqo. As antici- 
pated, this scaling yields good data collapse far from the 
edge, and even close to the edge. However, data does not 
collapse in the intermediate region, indicating a wave- 
dependent penetration length lp{q). 

Thus far, our estimates of the overall lengthscales of 
wavenumber amplification have shed no light on the na- 
ture of the cascade itself. The classic example of an elas- 
tic cascade was predicted by Pomeau and Rica Q for 
the so-called "curtain geometry" of a tension- free sheet 
rippled under a compressive force but constrained to be 
flat at one edge. They showed that the matching of the 
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FIG. 4: (A) A magnified image of the cascade. (B) At each 
value of X, a histogram of the scaled separation, god/(27r), 
between crests, for several values of distance x from the edge. 
Data were collected from two films with t = 246nm. The 
separations d, are determined from the locations of maxima 
of the intensity in the y direction. 

ripples to the flat edge (^e — > cx)) was achieved by an infi- 
nite hierarchy of branching events, in which each wrinkle 
branches into a succession of sharp folds with flat faces. 

One superficial difference between what we observe and 
the Pomeau-Rica cascade is that our cascade terminates 
at a finite wavenumber, and therefore passes through 
only a few generations. The more profound difference 
stems from the fact that our sheets experience a tension 
along the uncompressed direction. Any deviation from 
a one-dimensional wrinkling pattern imposes curvature 
in both directions; this Gaussian curvature generates in- 
plane stretching energy controlled by a modulus Y = Et. 
In the tensionless Pomeau-Rica scenario, the dominant 
contribution to the strain energy is the anharmonic en- 
ergy density ug ~ ^CxCyj whose minimisation leads to 
localised Gaussian curvature along a sequence of sharp 
ridges However, the consequence of the applied ten- 
sion is that the focusing of Gaussian curvature does not 
fully relieve the strain energy at other points: as noted 
in Eq. [31 the tension term ut penalises slope, and is 
non-zero even on flat facets where the gaussian curvature 
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vanishes. This mechanism thus favours a smooth reduc- 
tion of the ampUtude, which is naturally accomplished 
by the superposition of a finite number of Fourier modes 
with distinct wavenumbers. A transition between the 
Pomeau-Rica cascade and a smooth pattern is expected 
if ut ~ ug, which implies [see Supplementary Material] 
A < T/Y. Together with the threshold condition for 
wrinkling A > a^/Y we see that a necessary condition 
for a smooth pattern is that the tension be large enough 
such that £ ^ 1. The characteristic values of e in our ex- 
periments are well within this regime indicating a novel, 
smooth hierarchy, markedly different from the Pomeau- 
Rica stress-focusing cascade. 

A closer look at the cascade, as shown in the magnified 
view of Fig. 4A, supports this new scenario of a smooth 
mechanism in which larger amplitudes of higher wave- 
number Fourier components are smoothly mixed in as 
one approaches the edge. In Fig. 4B we present a more 
quantitative measure of the smoothness of the cascade. 
At a given distance x from the uncompressed edge we de- 
termine from the image, the separations d, between the 
crests of the wrinkles. At each value of x, we show a his- 
togram of Qod/ {2t:), the normalised separation between 
wrinkles. Far away from the edge, the separations are all 
concentrated at qod/{2TT) = 1. As expected, closer to the 
edge, more crests are formed, and at smaller values of 
d. Importantly, none of the histograms show significant 
weight near d = 0. In a scenario where wrinkles divide by 
localised branching, one might expect a preponderance of 
small values of d just after a branch-point, between sib- 
ling branches of the same parent wrinkle. That appears 
not to be the case in Fig. 4, with separations flowing 
smoothly to a mixture of higher Fourier components. 

The three forces operative in this problem - gravity, 
bending energies, and surface energies - can be combined 
in pairs to yield three distinct length-scales: {pg/By/^, 
(7/B)i/2 and Ic ^{l / pg)'^''^ ■ The first of these is qo, the 
second, an elastocapillary length [31 which controls qe, 
and the third is the capillary length, which determines 
the length of the cascade. However, we emphasise that 
the stress ratio e = cr/7 dictates the overall morphology 
of the pattern in our experiments. Surface tension plays 
a dual role in our experiments, determining both the en- 
ergy of the fluid meniscus, C/cap, as well as the tension ap- 
plied in the uncompressed direction. In principle, these 
are different effects that could be independently tuned. 
Increasing the capillary energy cost of the edge can tune 
a transition from the regime of our experiment to that 
of stiffer sheets, in which the effect of the edge is small. 
On the other hand, decreasing the applied tension could 
drive a transition from the smooth, energy-delocalised 
cascades we observe, to a regime of localised branching 
with energy- focusing. Thus, our observations open 



from solutions of polystyrene (PS; atactic, number- 
average molecular weight M„ = 121K, weight-average 
molecular weight = 1.05M„, radius of gyration 

Rg lOnm) in toluene, spin-coated on to glass substrates. 
The film thickness t was varied by changing the concen- 
tration of the solution and the spin rate, and was mea- 
sured by x-ray reflectivity using a Panalytical X-Pert x- 
ray diffractometer. The CuKa radiation from the x- 
ray source (wavelength A ~ 0.154nm) is coupled to a 
parabolic, graded multilayer mirror assembly that pro- 
duces a low- divergence beam of x-rays incident on the 
PS film. Observation of the Kiessig fringes in specular 
reflection with a two-circle goniometer yielded the film t 
with a precision of ±0.5nm. 

A rectangle of dimensions L x W — 2cm x 3cm was 
scribed onto the film with a sharp edge. When the sub- 
strate was dipped into a petri dish of distilled, deionized 
water, a rectangular piece of the PS film detached from 
the substrate. Because PS is hydrophobic, the film re- 
mained at the air-water interface. In order to prepare 
extremely well-defined edges, in a few selected cases, we 
spin-coated the film on a silicon wafer with an oxidised 
surface layer. Rather than scribing the film, the silicon 
substrate was fractured along a crystalline plane, and the 
PS film was lifted off with HF solution. 
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the way to the exploration of a rich phase diagram 
of both branched and smooth structures. 

Materials and Methods: The films were prepared 
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